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The DGLAP improved saturation model of Bartels, Golec and Kowalski is supplemented by 
the contribution from the heavy quarks: charm and beauty. It is shown to give good descrip- 
tion of both the total proton structure function F2 and, as a prediction, the heavy quarks 
contribution F%° and _F| 6 . The reasonable agreement for the longitudinal and diffractive 
structure functions is also found. 



Among a number of the dipole models of deep inelastic scattering the one proposed by 
Golec-Biernat and WusthofI™™ turned out to be especially popular. This is because it seems 
to grasp the essential elements of parton saturation being at the same time relatively simple. 
It gave a successful fit to low x HERA data for the structure function F%. However, with the 
advent of new more precise data the model required some improvement. This has been done by 
Bartels, Golec-Biernat and Kowalski^ (BGK model) by taking into account DGLAP evolution 
of the gluon density in the proton, which affects the behavior of the dipole cross section for small 
dipole sizes r. 

The photon-proton interaction at low x my be divided into two stages. In the first stage the 
transversely or longitudinally polarized virtual photon splits into a quark-anti-quark pair (a color 
dipole). This splitting is described by the photon wave function \P(r, z, Q 2 ,m 2 ), where z denotes 
the fraction of the photon momentum carried by quark in the dipole. In the second stage the 
dipole interacts with the proton with the dipole cross section a{x,r) which is usually modeled 
using ideas from high density QCD^l. In such description F2 is a sum of the contributions i 7 / 
coming from different quark flavours in the color dipole 

F 2 (x,Q 2 )=£f/(x,Q 2 ) (1) 
/ 



and 



Hi.x,Q 2 ) = -r-T— £ ld 2 rf l dz \¥ P (r, z,Q 2 ,m})\ 2 a (x,r) (2) 
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where P is the virtual photon polarization. 

In the BGK model the dipole cross section has the form 



a(x,r) =a l 1 - exp l-^—r 2 a s (n 2 )xg(x,fi 2 ) ]} (■'>) 



Table 1: Parameters of the DGLAP improved saturation model with heavy quarks obtained from the fit to F%. 
The results from the BGK fit with light quarks only are shown for reference in the first row. 
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with the scale \j? assumed to run with r 2 as /x 2 = C/r 2 + /xq and the gluon density xg(x,n 2 
which evolves with /j 2 according to the DGLAP equations with the initial condition 



xg(x,Ql) 
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at Ql = 1 GeV 2 



(4) 



The model has five parameters a^, C, /ig, A g and X g , which are fitted to the F2 data. In the 
paper of Bartels Golec-Biernat and Kowalski^ only the three light quarks were considered in 
the sum (^Q) . Thus the main goal of the work presented here is to consider the charm and beauty 
contributions within the BGK model. Notice that after the model parameters have been fitted 
to the F2 data, the contributions F$ 5 and F^ will come out as a prediction. 

One should emphasize that the heavy quarks production in deep inelastic scattering at low x 
could by no means be neglected. According to Hl^EJlTJ and zeusEI the contribution of the 
charm structure function to the total structure function varies from 10% to 30% while for the 
beauty structure function it reaches up to 3%. 

We performed fits with heavy quarks using the most recent HlEland ZEUS!!™ 

inclusive 

DIS data with x B < 0.01 and 0.1 < Q 2 < 4000 GeV 2 . We took the light quarks massless and the 
typical values of masses for heavy quarks: m c h arm = 1.3 GeV, mb ea uty = 5.0 GeV and modified 
the Bjorken variable to be x = x_b(1 + 4m 2 /Q 2 ), more details are given elsewhere^! The fit 
results are presented in Tabled In the first row we recall the parameters obtained in the BGK 
fit^ with light quarks only. Our fits with heavy quarks are presented in the next two lines. 
As we see from the value of x 2 /n&i the overall fit quality is still good. However, the values of 
parameters concerning the gluon differ significantly from the light quarks fit. In particular, X g 
is positive which means that the initial gluon distribution @ grows with decreasing x, which is 
opposite to the case with light quarks only when this distribution is valence-like. 

With the obtained parameters we are able to give predictions for the charm and beauty 
structure functions F^ c and F^ '. The results are plotted in Figs. ^ and |2 (solid lines). As a 
reference we show also the old GBW results without DGLAP modifications (dashed lines). We 
observe good agreement both in the normalization and in the slope for all Q 2 bins. The latter 
is improved by the DGLAP evolution, which is especially important at high Q 2 . 

The model provides also the predictions for the longitudinal structure function Fl and 



the diffractive structure function F, 



D(3) 



The obtained values of Fl are in agreement with HI 



estimations^ EOHl The comparison of this predictions with direct measurements, which are 
planned by HI and ZEUS collaborations, will be particularly interesting. For the case of F^ 
we still obtain reasonable agreement with the ZEUS data^l 

In summary, we have shown that the DGLAP improved saturation model with heavy quarks 
gives successful description of F2 at low x and provides predictions for the charm and beauty 
structure fun ctions F$ d and F$ b , which nicely agree with the data. For more details we refer to 
the paperQ2I. 
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Figure 1: Predictions for the charm structure function _F| C in the DGLAP improved saturation model (solid lines). 
The predictions in the old GBW model without DGLAP evolution are shown for reference (dashed lines). 




Figure 2: The same as in Fig.^but for the beauty structure function Ftf . 



